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Possible detection of relic neutrinos and their mass 

A. Ringwald 

Deutsches Elektronen- Synchrotron DESY, Notkestr. 85, D-22607, Hamburg, Germany 



Abstract. Recently the possibility was widely discussed that 
a large fraction of the highest energy cosmic rays may be 
decay products of Z bosons which were produced in the res- 
onant annihilation of ultrahigh energy cosmic neutrinos on 
cosmological relic neutrinos. If one takes this so-called Z- 
burst scenario seriously, one may infer the mass of the heav- 
iest relic neutrino as well as the necessary ultrahigh energy 
cosmic neutrino flux from a comparison of the predicted Z- 
burst spectrum with the observed cosmic ray spectrum. 



1 Introduction 

Big bang cosmology predicts the existence of a background 
gas of free photons and neutrinos. The measured cosmic mi- 
crowave background radiation (CMBR) supports the appli- 
cability of standard cosmology back to photon decoupling 
which occured approximately one hundred thousand years 
after the big bang. The relic neutrinos, on the other hand, 
have decoupled when the universe had a temperature of 1 
MeV and an age of just one second. Thus, a measurement of 
the rehc neutrinos, with a predicted average number density 
of 



56 cm 



(1) 



CMBR 



per light neutrino species i (m^. < 1 MeV), would provide 
a new window to the early universe. Their predicted number 
density is comparable to the one of the microwave photons. 
However, since neutrinos interact only weakly, the relic neu- 
trinos have not been detected untU now. 

Recently, an indirect detection possibility for relic neutri- 
nos has been discussed (Fargion et al., 1999; Weiler, 1999). 
It is based on so-called Z-bursts resulting from the resonant 
annihilation of ultrahigh energy cosmic neutrinos (UHECi^s) 
with relic neutrinos into Z bosons (Weiler, 1982; Roulet, 

Correspondence to: A. Ringwald (ringwald@mail.desy.de) 



1993; Yoshida et al. , 1997), the force carriers of the electro- 
weak neutral current (cf. Fig. 1). On resonance, the corre- 
sponding cross section is enhanced by several orders of mag- 
nitudes. If neutrinos have non-vanishing masses m^. - for 
which there is rather convincing evidence in view of the ob- 
servations of neutrino oscillations (Groom et al., 2000) - the 
resonance energies, in the rest system of the relic neutrinos, 
correspond to 



El: 



Ml 
2 m^^ 



4.2 • 10^1 eV 



1 eV 



(2) 



with Mz denoting the mass of the Z boson. These reso- 
nance energies are, for neutrino masses of 0(1) eV, remark- 
ably close to the energies of the highest energy cosmic rays. 
Indeed, it was argued (Fargion et al., 1999; Weiler, 1999) 
that the ultrahigh energy cosmic rays (UHECRs) above the 
predicted Greisen-Zatsepin-Kuzmin (GZK) cutoff (Greisen, 
1966; Zatsepin and Kuzmin, 1966; Bhattacharjee and Sigl, 
2000) around 4 • 10^^ eV are mainly protons (and, maybe, 
photons) from Z decay. 
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Fig. 1. Illustration of a Z-burst resulting from the resonant annihila- 
tion of an ultrahigh energy cosmic neutrino on a relic (anti-)neutrino 
(adapted from Ref. (Pas and Weiler, 2001)). 
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This hypothesis was discussed in several papers (Waxman, 
1998; Yoshida et al., 1998; Gelmini and Kusenko , 1999, 
2000; Blanco-Pillado et al., 2000; Pas and Weiler, 2001; Far- 
gion et al., 2001; McKellar et al. , 2001). Here, we review 
the Z-burst scenario and report on a recent quantitative in- 
vestigation (Fodor et al., 2001a, in prep.), where an attempt 
was made to determine the mass of the heaviest relic neutrino 
as well as the necessary UHECi^ flux by comparing the pre- 
dicted cosmic ray spectrum from Z-bursts with the observed 
one via a maximum likelihood analysis. 



2 Z-burst spectrum 

A comparison of the Z-burst scenario with the observed ul- 
trahigh energy cosmic ray spectrum proceeds essentially in 
four steps. First, one has to determine the probability of Z 
production as a function of the distance from Earth. Sec- 
ondly, one exploits data from collider experiments to derive 
the energy distribution of the produced protons in the lab sys- 
tem. Thirdly, one has to take into account the propagation of 
the protons to Earth, i. e. one has to determine their energy 
losses due to pion and e+e~ production through scattering 
on the CMBR and due to their redshift. The last step is the 
comparison of the predicted and observed spectra and the ex- 
traction of the mass of the relic neutrino and the necessary 
VHECv flux. 

The contribution of protons from Z-bursts to the UHECR 
spectrum can be summarized as (Fodor et al., 2001a, in prep.) 

oo Ro oo 

* 

X [F,, {E,^ , r) np. (r) + Fp, {E^, , r) n,, (r)] 
XCT^iPi(e) Bt{Z hadrons) Qp+n{Ep) 

x{^)—Pir,Ep;E), 

where E is the energy of the protons arriving at Earth. Fur- 
ther important ingredients in Eq. (3) are: the ultrahigh energy 
cosmic neutrino fluxes F,^. [E^. , r) at the resonant energies 
E^i ~ E]J^^ and at distance r to Earth, the number density 
rii,. ir) of the relic neutrinos, the Z production cross section 
CTi^.p. (e) at centre-of-mass (cm) energy e = 2 rriu^ E^. , the 
branching ratio for Z to hadrons, Bi{Z hadrons), the 
energy distribution Qp+n{Ep) of the produced protons (and 
neutrons) with energy Ep, and the probability P{r, Ep; E) 
that a proton created at a distance r with energy Ep arrives at 
Earth above the threshold energy E. 

The last four building blocks, CTiy^p;, the hadronic branch- 
ing ratio, Q, and P, are very well determined, whereas the 
first two ingredients, the flux of UHECi^s, F^. [Ei,. , r), and 
the relic neutrino number density n^. (r), are much less ac- 
curately known. In the following we shall discuss all these 
ingredients in detail. 
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Fig. 2. Proton plus antiproton momentum distribution in Z de- 
cays (from Ref. (Fodor et al., in prep.)). Top: Combined data from 
Ref. (Akers et al., 1994; Abreu et al., 1995; Buskulic et al., 1995; 
Abe et al., 1999, unpublished). Bottom: Distribution in the lab sys- 
tem, in which the target relic neutrino is at rest. 



2.1 Z production and decay 

At LEP and SLC millions of Z bosons were produced and 
their decays analyzed with extreme high accuracy. Due to 
the large statistics, the uncertainties of the analysis related to 
Z decay are negligible. 

Fodor et al. (2001a, in prep.) combined existing published 
(Akers et al., 1994; Abreu et al., 1995; Buskulic et al., 1995; 
Abe et al., 1999) and some improved unpublished (Abe et al., 
unpublished) data on the momentum distribution Vp (x) = 
dNp/dx, with momentum fraction x = Pproton/pbcam, of 
protons (p) (plus antiprotons (p)) in Z decays, see Fig. 2 
(top). The p + p multipHcity is (Np) = dxVp{x) = 
1.04 ± 0.04 in die hadronic channel (Groom et al., 2000). 

The energy distribution Qp{Ep) of the produced protons 
with energy Ep entering in the spectrum (3) is obtained after 
a Lorentz transformation from the cm system to the lab sys- 
tem, in which the target relic neutrino is at rest. It is only a 
function of y — 2Ep/Ei, and displayed in Fig. 2 (bottom). 

Neutrons produced in Z decays will decay and end up as 
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Fig. 3. Upper limits on differential neutrino fluxes in the ul- 
trahigh energy regime. Shown are experimental upper limits on 
Fjj^ + Fv^ from Fly's Eye (Baltrusaitis et al., 1985) and on 
X^j^^g J\Fvf + Pvf) from the Goldstone lunar ultrahigh energy 
neutrino experiment GLUE (Gorham et al., 2001), as well as theo- 
retical upper limits on F^^ + F^^ from "visible" ("WB" (Waxman 
and Bahcall , 1999; Bahcall and Waxman , 2001)) and "hidden" 
("MPR" (Mannheim et al. , 2001)) hadronic astrophysical sources. 



UHECR protons. They are taken into account according to 

(4) 



where the neutron (n) (h- antineutron (n)) multiplicity, (N^), 
is w 4% smaller than the proton's (Biebel , priv. comm.). 

2.2 Propagation of nucleons through the CMBR 

Similarly, the CMBR is known to a high accuracy. It plays 
an important role in the determination (Bahcall and Waxman, 
2000; Fodor and Katz, 2001a) of the probabiHty P(r, E^,, E), 
which takes into account the fact that protons of extragalactic 
(EG) origin and energies above « 4-10^^ eV lose a large frac- 
tion of their energies (Greisen, 1966; Zatsepin and Kuzmin, 
1966) due to pion and e+e^ production through scattering 
on the CMBR and due to their redshift. P(r, Ep,E), in the 
form as it has been calculated for a wide range of param- 
eters by Fodor and Katz (2001a), is an indispensible tool 
in the quantitative analysis in Refs. (Fodor et al., 2001a, in 
prep.). In addition to its already published (Fodor and Katz, 
2001a) parametrized form, the numerical data for the prob- 
ability distribution {—)dP{r,Ep;E)/dE are now available 
via the World- Wide- Web URL www.desy.de/~uhecr . 

2.3 UHECi^ fluxes 

Presently unkown ingredients in the evaluation of the Z-burst 
spectrum (3) are the differential fluxes F^. of ultrahigh en- 
ergy cosmic neutrinos (see e.g. Refs. (Protheroe , 1999; 
Gandhi , 2000; Learned and Mannheim , 2000) for recent 
reviews). Present experimental upper limits on these fluxes 
are rather poor (cf. Fig. 3 and contributions to these proceed- 
ings (Seckel et al. , 2001; Yoshida et al. , 2001)). 



What are the theoretical expectations? More or less guar- 
anteed are the so-called cosmogenic neutrinos which are pro- 
duced when ultrahigh energy cosmic protons scatter inelasti- 
cally off the cosmic microwave background radiation (Grei- 
sen, 1966; Zatsepin and Kuzmin, 1966) in processes such as 
P7 — > A — ^ n7r+, where the produced pions subsequently 
decay (Berezinsky and Zatsepin , 1969, 1970). These fluxes 
(for recent estimates, see (Yoshida and Teshima , 1993; Pro- 
theroe and Johnson , 1996; Yoshida et al. , 1997; Engel and 
Stanev , 2001)) represent reasonable lower limits on the ul- 
trahigh energy cosmic neutrino flux, but turn out to be in- 
sufficient for the Z-burst scenario. Recently, theoretical up- 
per limits on the ultrahigh energy cosmic neutrino flux have 
been given by Waxman and Bahcall (1999), Mannheim et 
al. (2001), and Bahcall and Waxman (2001). Per con- 
struction, the upper limit from "visible" hadronic astrophys- 
ical sources, i. e. from those sources which are transparent 
to ultrahigh energy cosmic protons and neutrons, is of the 
order of the cosmogenic neutrino flux (Waxman and Bah- 
call , 1999; Mannheim et al. , 2001; Bahcall and Waxman , 
2001) and shown in Fig. 3 ("WB"). Also shown in this fig- 
ure ("MPR") is the much larger upper limit from "hidden" 
hadronic astrophysical sources, i. e. from those sources from 
which only photons and neutrinos can escape (Berezinsky , 
1979; Mannheim et al. , 2001). 

In this situation of insufficient knowledge, the following 
approach concerning the flux of ultrahigh energy cosmic neu- 
trinos was taken in Refs. (Fodor et al., 2001a, in prep.). The 
flux was assumed to have the form 



(5) 



where z is the redshift and where a characterizes the source 
evolution (see also (Yoshida et al. , 1997; Yoshida et al., 
1998)). The flux at Earth, F,j.{E,y.,0), or, more accurately, 
the sum of the corresponding fluxes at the resonant energies. 



Fz 



E 



F,,{El^^) + F.^iEl^^)] 



(6) 



was then determined, together with the neutrino mass, by a 
fit to the UHECR data. This analysis went up to distances 
Rq (cf. (3)) corresponding to redshift z = 2 (cf. (Waxman, 
1995)), and uncertainties in the Hubble expansion rate Hq = 
100 h kms^i Mpc"\ as given in Ref. (Groom et al., 2000), 
were included. 

2.4 Relic neutrino number density 

The dependence of the relic neutrino number density n,y. on 
the distance r was treated in the following way in Refs. (Fo- 
dor et al., 2001a, in prep.). The question is whether there is 
remarkable clustering of the relic neutrinos within the local 
GZK zone of about 50 Mpc. It is known that the density 
distribution of relic neutrinos as hot dark matter follows the 
total mass distribution; however, with less clustering (Ma , 
1999; Primack and Gross , 2000). In fact, for <1 eV, 
one expects pretty much that the neutrino number density 
equals the big bang prediction (1) (Primack , priv. comm.). 
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Fig. 4. Mass density fluctuation field 5 along the supergalactic plane 
as obtained from peculiar velocity measurements (da Costa et al. , 
1996). Shown are contours in intervals of 5 = 2, surface maps 
on a grid of spacing 500 kms~^, corresponding to 5h~^ Mpc, 
with the height proportional to 5, and contrast maps. One recog- 
nizes some well-known structures in the nearby volume such as 
the Great Attractor at supergalactic coordinates (SGX ~ —2000 
kms~^, SGY ~ —500 kms~^), the Perseus-Pisces complex (SGX 
~ 6000 kms"\ SGY ~ -1000 kms"^), and the large void (SGX 
~ 2500 kms"\ SGY ~ kms"^) in between. One can also see 
traces of the Coma/A 1 367 supercluster along SGY ^ 7000 km s~ ^ , 
and the Cetus region (SGX ~ 500 km s"\ SGY ~ -6000 kms'^). 



To take above facts into account, the shape of the (r) dis- 
tribution was varied in Refs. (Fodor et al., 2001a, in prep.), 
for distances below 100 Mpc, between the standard cosmo- 
logical homogeneous case (1) and that of mtot(f), the total 
mass distribution obtained from peculiar velocity measure- 
ments (da Costa et al. , 1996; Dekel et al., 1999) (cf. Fig. 4). 
These peculiar measurements suggest relative overdensities 
of at most a factor 2-^3, depending on the grid spacing (cf. 
Fig. 4). A relative overdensity of 10^ ~ 10^ in our neigh- 
bourhood, as it was assumed in earlier investigations of the 
Z-burst hypothesis (Fargion et al., 1999; Weiler, 1999; Wax- 
man, 1998; Yoshida et al., 1998; Blanco-Pillado et al., 2000) 
(see also Ref. (McKellar et al. , 2001)), seems unnatural in 
view of these data. The quantitative results turned out to be 
rather insensitive to the variations of the overdensities within 
the considered range, whose effect were included in the final 
error bars. For scales larger than 100 Mpc the relic neutrino 
density was taken according to the big-bang cosmology pre- 
diction, n^. = 56 ■ (1 + z)'^ cm^"^. 



3 Determination of and the UHECi^ flux 

The predicted spectrum of protons from Z-bursts (3) can now 
be compared with the observed ultrahigh energy cosmic ray 
spectrum (cf. Fig. 5). The analysis in Refs. (Fodor et al., 
2001a, in prep.) included published and unpublished (from 
the World-Wide- Web pages of the experiments on 17/03/01) 
UHECR data of AGASA (Takeda et al., 1998, 1999), Fly's 
Eye (Bird et al., 1993, 1994, 1995), Haverah Park (Lawrence 
et al., 1991; Ave et al., 2000), and HIRES (Kieda et al., 
2000). Due to normalization difficulties the Yakutsk (Efimov 
et al., 1991) results were not used. 

As usual, each logarithmic unit between log(_E/eV) = 18 
and log(_E/eV) — 26 was divided into ten bins. The pre- 
dicted number of ultrahigh energy cosmic ray events in a bin 
was taken as 

N{t) = (7) 
dE [f;,|bkd(S; A, 13) + Fp\z{E; Fz,m,.)\ , 



Ei 



where Ei is the lower bound of the i energy bin. The 
first, background term Fp|bkd(-E; A, /3) was taken to have 
the usual power-law behavior which describes the data well 
for smaller energies (Takeda et al., 1998, 1999). The second 
term of the flux in Eq. (7) corresponds to the spectrum of the 
Z-bursts, Eq. (3). 

Two possibilities for the background term were studied. 
In the first case it was assumed that the power part is pro- 
duced in our galaxy. Thus no GZK effect was included for it 
("halo"), and it was taken as. 



Fp\^i,^{E-A,P)=A-E 



-0 



(halo) , 



(8) 



In the second - in some sense more realistic - case it was 
assumed that the background protons come from uniformly 
distributed, extragalactic sources and suffer from the GZK 
cutoff ("EG"). In this case, A ■ E^ ^ was taken as an injection 
spectrum, and the simple power-law-like term was modified, 
by taking into account the probability P{r, Ep \ E), 



Fp\hkdiE; A, /3) 



X A ■ E^P (-) 



dEr, 



dr 



(9) 



dPir,Ep;E) 
dE 



(EG) 



It falls off around 4 • IQ^*^ eV (see Fig. 5 (bottom)). 

Note that the following implicit assumptions have been 
made through ansatz (7): ;) It was assumed that the ultra- 
high energy photons from Z-bursts can be neglected. //) It 
was assumed that there are no significant additional primary 
ultrahigh energy proton fluxes beyond the extrapolation of 
the above power-law. Assumption i) can be justified on ac- 
count of the result of a detailed study (Fodor et al., in prep.) 
of the boosted Z-decay (data from Refs. (Akers et al., 1994; 
Abreu et al., 1995; BuskuUc et al., 1995; Abe et al., 1999, 
unpublished)) which show that the energy of the photons is 
peaked at 1.7 • 10^^ (1 eNIm^.) eV, at which the attenuation 
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Fig. 5. The available UHECR data with their error bars and the best 
fits from Z-bursts, for nearly degenerate neutrino masses, m^^ « 
and for cosmological evolution parameter a = 3 (Fodor et 
al., 2001a, in prep.). Top: Best fit for the "halo"-case (solid line). 
The bump around 4 ■ 10^^ eV is due to the Z-burst protons (dash- 
dotted), whereas the almost horizontal contribution (dashed) is the 
first, power-law-like term of Eq. (7). Bottom: "Extragalactic"-case. 
The first bump at 4 ■ lO"' eV represents protons produced at high 
energies and accumulated just above the GZK cutoff due to their 
energy losses. The bump at 3 • 10^^ eV is a remnant of the Z-burst 
energy. The dashed line shows the contribution of the power-law- 
like spectrum with the GZK effect included. The predicted fall-off 
for this term around 4 ■ lO"' eV can be observed. 

length of photons is much smaller than that of protons. Thus, 
the contribution of ultrahigh energy photons from Z decay in 
the observed UHECR spectrum is far less relevant than that 
of the protons. Assumption ii) will be relaxed in Ref. (Fodor 
et al., in prep.). 

The expectation value for the number of events in a bin 
is given by Eq. (7). To determine the most probable value 
for nil, the maximum likelihood method was used and the 
X^(/3, A, Fz, m^) (Fodor and Katz, 2001b), 

26.0 

= ^ 2 [N{i} - Nod) + iVo(i) In (iVo(*)/iV(z))] ,(10) 

1=18.5 

was minimized, where No{i) is the total number of observed 
events in the i'^ bin. Since the Z-burst scenario results in 
a quite small flux for lower energies, the "ankle" was used 
as a lower end for the UHECR spectrum: log(i?niin/eV) = 



18.5. The results are insensitive to the definition of the up- 
per end ( - the flux is extremely small there - ) for which 
log(i?„iax/eV) = 26 was chosen. The uncertainties of the 
measured energies are about 30% which is one bin. Using 
a Monte-Carlo method this uncertainty was included in the 
final error estimates. 

In the case where there are three neutrino types with nearly 
degenerate neutrino masses, m^. « nii,, the fitting procedure 
involves four parameters: A, Fz and m^. The minimum 
of the x^(/3, A, Fz, m^) function is x^^^, at m^^-ain which is 
the most probable value for the mass, whereas the 1 a (68%) 
confidence interval for is determined by 

X\P', A', F'z,m,) = xl{m,) = xLn + 1 ■ dD 

Here P', A', F^ are defined in such a way that the function 
is minimized in /3, A and Fz, at fixed m^. 

Qualitatively, the analysis can be understood in the follow- 
ing way. In the Z-burst scenario a small relic neutrino mass 
needs large El°^ in order to produce a Z. Large El°^ results 
in a large Lorentz boost, thus large Ep. In this way the shape 
of the detected energy (E) spectrum determines the mass of 
the relic neutrino. The sum of the necessary UHECi^ fluxes 
Fz, on the other hand, is determined by the over-all normal- 
ization. 

The best fits to the observed data can be seen in Fig. 5, 
for evolution parameter a = 3. For the "halo"-case, a neu- 
trino mass of 2.75^Qgy|^ggj eV was found, whereas, for the 

"EG"-case, the fit yielded 0.26^Q^|i||[J'22j eV. The first num- 
bers are the 1 a, the numbers in the brackets are the 2 a er- 
rors. This gives an absolute lower bound on the mass of 
the heaviest v of 0.06 eV at the 95% confidence level (CL). 
Note, that the surprisingly small uncertainties are based on 
the above x^ analysis and dominantly statistical ones. The 
fits are rather good; for 21 non-vanishing bins and 4 fitted 
parameters they can be as low as x^ — 18.6. The neu- 
trino mass was determined for a wide range of cosmolog- 
ical source evolution (a = 0-^3) and Hubble parameter 
(TJq = (71 ± km/sec/Mpc) and was found to de- 

pend only moderately on them. The results remained within 
the above error bars. A Monte-Carlo analysis was performed 
studying higher statistics. In the near future. Auger (Guer- 
ard, 1999; Bertou et al., 2000) will provide a ten times higher 
statistics, which reduces the error bars in the neutrino mass 
to K, one third of their present values. 

The necessary UHECj^ fluxes at E]^^ have been obtained 
from the fits. They are summarized in Fig. 6, together with 
some existing upper limits and projected sensitivities of pre- 
sent, near future and future observational projects. The nec- 
essary UHECi^ flux appears to be well below present upper 
limits and is within the expected sensitivity of AMANDA, 
Auger, and OWL. An important constraint for all top-down 
scenarios (Bhattacharjee and Sigl, 2000) is the EGRET ob- 
servation of a diffuse 7 background (Sreekumaretal. , 1998). 
As a cross check, one may calculate the total energy in pho- 
tons from Z-bursts. Even if one assumes that all their energy 
ends up between 30 MeV and 100 GeV, one finds that the 7 



6 



r 10'° r ^ GLUE 

- 10" r Fiys Eye 



El. 

w 
H 



I I I iiiii| — I I I iiiii| — I I I iiiii| — I I I iiiii| — I I I iiiiii — I I 1 1 




"halo'.v"^^' 



^q4 t- OWL 

' I I I I I ^ 




10" 10'" 10'° 10" 10'' 
E [eV] 

Fig. 6. Neutrino flux, -F = | Z^Li (-^"^^ + Fi^,)' required by the Z- 
burst hypothesis for the "halo" and the "extragalactic" case, for evo- 
lution parameter a = 3 (Fodor et al., 2001a, in prep.). The horizon- 
tal errors indicate the 1 a uncertainty of the mass determination and 
the vertical errors include also the uncertainty of the Hubble expan- 
sion rate. The dependence on a is just of the order of the thickness 
of the lines. Also shown are upper limits from Fly's Eye (Baltru- 
saitis et al., 1985) and the Goldstone lunar ultrahigh energy neutrino 
experiment GLUE (Gorham et al., 2001), as well as projected sen- 
sitivities of AMANDA (Barwick , unpublished). Auger (Yoshida 
et al., 1998; Capelle et al., 1998) and OWL (Yoshida et al., 1998; 
Ormes et al., 1997). 



flux is somewhat smaller than that of EGRET. These numer- 
ical findings are in fairly good agreement with simulations 
done by Sigl (priv. comm.) (cf. Fig. 7 (top)). 

4 Comparison with Amf, from neutrino oscillations 

One of the most attractive patterns for neutrino masses is si- 
milar to the one of the charged leptons or quarks: the masses 
are hierarchical, thus the mass difference between the fami- 
lies is approximately the mass of the heavier particle. Using 
the mass difference of the atmospheric neutrino oscillation 
for the heaviest mass (Groom et al., 2000), one obtains val- 
ues between 0.03 and 0.09 eV. It is an intriguing feature of 
the result found in Refs. (Fodor et al., 2001a, in prep.) that 
the smaller one of the predicted masses is compatible on the 
w 1.3 (T level with this scenario. 

Another popular possibility is to have 4 neutrino types. 
Two of them - electron and sterile neutrinos - are separated 
by the solar neutrino oscillation solution, the other two - 
muon and tau - by the atmospheric neutrino oscillation solu- 
tion, whereas the mass difference between the two groups is 
of the order of 1 eV. This possibility was studied, too. On the 
mass scales and resolution considered, the electron and ster- 
ile neutrinos are practically degenerate with mass mi and the 
muon and tau neutrinos are also degenerate with mass m2- 
The best fit and the 1 cr region in the mi — m2 plane is shown 
in Fig. 8 for the "EG"-case for the background protons. The 
dependence of this result on the cosmological evolution and 
on the UHECi^ spectrum will be discussed in Ref. (Fodor et 
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Fig. 7. Spectra of nucleons, 7 rays, and neutrinos for three degen- 
erate neutrinos m,^^ = = 1 eV, for cosmological evolution 
parameter a = 3 and an input differential spectrum of ultrahigh 
energy cosmic neutrinos oc E^^. Top: Obtained with the standard 
cosmological value (1) for the relic neutrino number density (Sigl , 
priv. comm.). Bottom: Obtained with a local overdensity of 300 on 
a scale of 5 Mpc over the the standard cosmological value (1) for 
the relic neutrino number density (Yoshida et al., 1998). 



al., in prep.). Since this two-mass scenario has much less 
constraints the allowed region for the masses is larger than in 
the one-mass scenario. 



5 Comparison with other studies 

Numerical simulations of Z-burst cascades for rrii, ^ 1 eV, 
taking into account all known extragalactic propagation ef- 
fects, were performed by Yoshida et al. (1998). Based on 
case studies, local overdensities by factors ranging from 20^ 
10"^ over the standard cosmological relic neutrino number 
density on a scale of 5 Mpc were argued to be necessary in 
order to get a successful description of the UHECR events 
and rate above the GZK cutoff without violating lower en- 
ergy photon flux limits and without invoking inconceivable 
UHEC;/ fluxes. For such large overdensities, most of the 
UHECRs from Z-bursts originate nearby and their attenua- 
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Fig. 8. The best fit and the 1 cr (68% confidence level) region in 
a scenario with two non-degenerate neutrino masses (Fodor et al., 
2001a, in prep.). 



tion to the Earth can be neglected. In our case, with real- 
istic overdensities < 2 3 on scales < 100 Mpc, most of 
the UHECRs from Z-bursts originate from cosmological dis- 
tances and suffer therefore much more from GZK attenuation 
effects. Therefore, despite of the fact that by construction 
the overall rate of UHECRs from Z-bursts observed at Earth 
is the same in both investigations, the predicted spectra are 
quite different, as illustrated in Fig. 7. For example, in the 
case of a large neutrino overdensity, the Z-burst proton spec- 
trum has a single bump, reflecting that of Fig. 2 (bottom), at 
some fraction of the resonance energy (cf. Fig. 7 (bottom)), 
whereas it is seriously deformed in the case with a realistic 
relic neutrino number density, piling up near the GZK cutoff 
(cf. Figs. 5 and 7 (top)). Moreover, whereas, in the over- 
dense case, the contribution of Z-burst photons dominates, in 
the presently observed ultrahigh energy region, over the one 
of Z-burst protons (cf. Fig. 7 (bottom)), it can be neglected 
in the non-overdense case (cf. Fig. 7 (top)). 

Both variants of the Z-burst scenario, the one with and the 
one without remarkable local clustering of reUc neutrinos, 
have problems. In the former case one has to admit that no 
mechanism is known which leads to such a large overden- 
sity. Gravitational clustering is really too inefficient for neu- 
trinos with Txiy < 0(1) eV (Ma , 1999; Primack and Gross , 
2000; Primack , priv. cormn.). The case without remarkable 
local overdensity, on the other hand, requires a tremendous 
ultrahigh energy cosmic neutrino flux (cf. Fig. 6), which, 
however, may be mildered by a factor of 0(3) if one in- 
vokes, in extension of the standard cosmology, lepton asym- 
metries (Gelmini and Kusenko , 1999) at a level which is 
still allowed (Kneller et al. , 2001). The required neutrino 
flux is larger than the theoretical upper limit from "hidden" 
hadronic astrophysical sources, from which only photons and 
neutrinos can escape (cf. Fig. 3 ("MPR")). From an analysis 
of the assumptions behind this limit (Mannheim et al. , 2001) 
one finds that one has to invoke hidden astrophysical sources 
whose photons somehow don't show up in the diffuse 7 ray 



background measured by EGRET (Sreekumar et al. , 1998). 
One of the possibilities are sources whose surrounding is so 
dense that no ultrahigh energy photons can escape. It is an 
interesting question whether such conditions can be reaUzed 
in BL Lacertae objects, a class of active galactic nuclei re- 
cently discussed as possible sources of the highest energy 
cosmic rays (Tinyakov and Tkachev , 2001). Alternatively, 
one may invoke top-down scenarios (Bhattacharjee and Sigl, 
2000) for the sources of the highest energy cosmic neutrinos 
such as unstable superheavy relic particle decays (Gelmini 
and Kusenko , 2000). 



6 Conclusions 

I reviewed a recent comparison of the predicted spectrum 
from Z-bursts - resulting from the resonant annihilation of 
ultrahigh energy cosmic neutrinos with relic neutrinos - with 
the observed ultrahigh energy cosmic ray spectrum (Fodor 
et al., 2001a, in prep.). The mass of the heaviest rehc neu- 
trino turned out to be my = 2.75lo'97 halo- and 
0.26^0 ^4 eV for extragalactic-scenarios for the background 
protons. The second mass, with a lower bound of 0.06 eV on 
the 95% CL, is compatible with a hierarchical neutrino mass 
scenario with the largest mass suggested by the atmospheric 
neutrino osciUation. The above neutrino masses are in the 
range which can be explored by future laboratory experi- 
ments like the f3 decay endpoint spectrum and the neutrino- 
less (3(3 decay (Pas and Weiler, 2001; Vissani , 2001; Farzan 
et al., 2001; Osland and Vigdel , 2001; Czakon et al. , 2001) 
or by simultaneous observations of the bursts of neutrinos 
and gravitational waves emitted during a stellar collapse (Ar- 
naud et al. , 2001). They compare also favourably with the 
upper limit rriiy^ < 4.4 eV found recently from a global 
cosmological analysis involving also the recent CMBR mea- 
surements (Wang et al., 2001), with a future sensitivity down 
to <0.3 eV (Hu et al. , 1998). We analysed a possible two- 
mass scenario and gave the corresponding confidence level 
region. The necessary UHECz^ flux was found to be consis- 
tent with present upper limits and detectable in the near fu- 
ture. Astrophysical sources of these ultrahigh energy cosmic 
neutrinos should be hidden even in 7 rays. 
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